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Abstract 

In light of the Higgs discovery and the nonobservation of sparticles at the LHC, we revisit the 
SUSY induced top quark flavor changing decay into the Higgs boson. We perform a scan over the 
relevant SUSY parameter space by considering the constraints from the Higgs mass measurement, 
the LHC search for SUSY, the vacuum stability, the precision electro-weak observables as well as 
B — X,y. We have the following observations: (1) In the MSSM, the branching ratio of t > ch 
can only reach 3.0 x 1075, which is about one order smaller than previous results obtained before 
the advent of the LHC. Among the considered constraints, the Higgs mass and the LHC search for 
sparticles are found to play an important role in limiting the prediction. (2) In the singlet extension 
of the MSSM, since the squark sector is less constrained by the Higgs mass, the branching ratio 
of t — ch can reach the order of 107? in the allowed parameter space. (3) The chiral-conserving 
mixings 6,7, and ógg may have remanent effects on t > ch in heavy SUSY limit. In the MSSM 
with squarks and gluino above 3 TeV and meanwhile the CP-odd Higgs boson mass around 1 TeV, 


the branching ratio of t — ch can still reach the order of 107? under the constraints. 


PACS numbers: 


I. INTRODUCTION 


A scalar with mass around 125 GeV has been discovered at the LHC a "n According 
to the analysis of the ATLAS and CMS collaborations, the measured properties of this 
scalar, albeit with large experimental uncertainties, agree well with those of the Higgs boson 
in the Standard Model (SM), which means that it plays a role in the electroweak (EW) 
symmetry breaking and also in the mass generation for the fermions in the SM . Even 
so, due to the deficiencies of the SM itself in describing the symmetry breaking, it is well 
motivated to interpret this scalar in various frameworks of new physics. Obviously, in 
order to ambiguously decipher the nature of the scalar, it is mandatory to scrutinize both 
experimentally and theoretically the couplings of the scalar, including its self-interactions. 
In this direction, the couplings of the scalar with the yet known heaviest particle, top quark, 
are of fundamental importance since, as suggested by the LHC Higgs data, the htt coupling 
is strong, and meanwhile it is widely conjectured to be sensitive to new physics. In fact, great 
efforts have been paid recently to investigate the top-Higgs associated production processes 
like pp — tth y and pp — qth at the LHC to extract the size and sign of the htt 
Yukawa coupling, and also the top quark flavor changing decay t — ch to prob anomalous 
top-Higgs interaction [9], |10]. 

Among the new physics models, the supersymmetric theory (SUSY) is a promising one 
due to its capability to solve the hierarchy problem of the SM, unify the gauge coupling as 
well as provide a viable Dark Matter candidate] . In SUSY, a SM-like Higgs boson 
h around 125GeV usually implies third generation squarks at or heavier than 1TeV, and 
the preference of the heavy squarks is further corroborated by the absence of any signal in 
the search for SUSY at the LHC. If the SUSY scale is really high, which was focused on 
in many recent theoretical works{13) the only way to detect SUSY is through its possibly 
large remanent effects in EW processes. Such effects may exist in the Higgs process because 
the dominant_part of the Higgs couplings to squarks is proportional to soft SUSY breaking 
parameters|11], and consequently, the suppression induced by the squark propagators in 
SUSY radiative correction to the process may be compensated under certain conditions. 
This feature has been demonstrated in the SUSY correction to the hbb vertex(14) the Higgs 
pair production process at the Lucha. and also the Higgs rare decay h — rallo. Here 


we emphasize that the existence of the remanent effect in the asymptotic large SUSY mass 


limit does not contradict with the Appelquist-Carazzone theorem| 17, which is valid only for 
supersymmetric theories with an exact gauge symmetry. Previous studies on such remanent 
effects in SUSY, also see [18]. 

In this work, we focus on the top quark flavor changing decay t — ch in SUSY. The 
reasons that we are interested in it mainly come from three considerations. Firstly, the LHC 
as a top factory has great capability to scrutinize the properties of top quark, including its 


rare decay modes. As far as the flavor changing decay t — ch is concerned, its branching 


ratio in the SM is only at the order of 107 !^|19], while in SUSY it may be greatly enhanced 
to 1074 according to previous studies " 21]. Since any observation of the decay in future 
will be a robust evidence of new physics, this decay should be paid attention to in the LHC 
era, especially noting the fact that the Higgs boson has been recently discovered. Secondly, 
as introduced before, the LHC experiment has measured Higgs mass and pushed SUSY to a 
rather high scale. These results have great impacts on the SUSY prediction about t — ch, 
so it is necessary to update previous studies on t > ch in light of the experimental progress. 
Thirdly, unlike the other top FCNC processes in uL the decay t — ch may have 
remanent effect in heavy SUSY limit. From theoretical point of view, it is worthwhile to 
investigate such a feature in detail. Besides, we remind that, if the flavor mixings between 
scharm and stop are present, which may push up the rate of t — ch greatly, the LHC 
constraint on stop masses can be relaxed. This in return can alleviate the fine tuning 
problem of the SUSY 23 : 

This paper is organized as follows. In section II, we parameterize the flavor mixings 
in squark sector and define our conventions. We also list various constraints on SUSY. In 
section III, we study the decay t — ch in both low energy SUSY and heavy SUSY, and 
present some benchmark points at which the predictions on t — ch are optimized. We also 
exhibit the features of the remanent effect on t > ch. Finally, we present our conclusions in 


section IV. 


II. FCNC INTERACTIONS IN SUSY 


In the supersymmetric theories such as the Minimal Supersymmetric Standard Model 
(MSSM)[LI] and the Next-to Minimal Supersymmetric Standard Model (NMSSM) Lig. the 


squark sector consists of six up-type squarks (ür, Čz, tL, üg, Čr, tg) and six down-type 


squarks (d L, ŠL, bL, dR, ŠR, br). In general, the states with different chiral and flavor quan- 
tum numbers in each type of squarks will mix to form mass eigenstates, and consequently 
potentially large flavor changing interactions arise from the misalignment between the rota- 


tions that diagonalize quark and squark sectors. In the super-CKM basis, the 6 x 6 squark 


mass matrix M? (q = à, d) takes the form |24 


(M2iLc-C;^ (Mĉ)er-— CF” (1) 
(QMj)ua — CEP)! (M2) nn + OF" 


à 
= 


where C^ = m? + cos26M3(T3 — Q,82,)1, CPR = m2 cos 28M2Q,sljl and CPR = 
mapu(tan B)-?13 are 3 x 3 diagonal matrices with 1 standing for the unit matrix in flavor 
space, m, being the diagonal quark mass matrix and T? = 1, -i for q = u, d respectively, and 
tan B = 22 is the ratio of the vacuum expectation values of the SU(2) doublet Higgs fields. If 
one only considers the flavor mixings between the second and the third generation squarks, 


the soft breaking squared masses (M?)1r, (Mz)rg and (M?)nn can be parameterized as 


Mà 0 0 
(Majte = 0 Mọ, rr Mo, Mg; |> 
0 ðLLMQ, MQ; Mé, 
0 0 0 
(MZ)rr = | 0 0 órgvoMpg | > 
0 Ónr US Mgr, ma A, 
(Mi)ng = (Mi)Lrlug, 5M$,. iion (2) 


where Mg, and My, (i = 1,2,3 denotes generation index) are soft breaking parameters with 
mass dimension, MP, and Mf, represents SUSY scale defined as Mpa = (Mu, + Mo;)/2 
and MẸ, = (My, + Mg,)/2, and órr, órg, Or, and Ogg reflect the extent of the flavor 
violation. Similarly, for down-squarks we have 
0 0 0 
(Miis = | 0 0 o gviMER |> 
0 ôt vM, MA 
(Mī)rr = (Mi)LLIu, >m., 114649) (3) 
and due to SU (2) gauge symmetry, (M3) LL is determined wb 
(Mj)ur = Végu(M$) 1L Vou (4) 


4 


with Vexm denoting the Cabibbo-Kobayashi-Maskawa matrix in the SM. Note that in 
Eqs.(2]3), we only keep the chiral-flipping terms for third-family squarks because these 
terms are usually assumed to be proportional to corresponding quark masses, and can not 
be neglected only for third-family squarks. 

The squark mass eigenstates can be obtained by diagonalizing the mass matrix presented 
above with an unitary rotation U;, which is performed numerically in our analysis. The 


interaction of the field X with a pair of squark mass eigenstates is then obtained by 
V(X4jy) = Ui Uje V(XEG) | (5) 


where V(X qq) denotes a generic vertex in the interaction basis and V (X q;, 


405) is the vertex 


in the mass-eigenstate basis. It is clear that both the squark masses and their interactions 
depend on the mixing parameters ó;s. 

In some fundamental supersymmetric theories like the mSUGRA and gauge-mediated 
SUSY-breaking models, the mixing parameters are functions of the soft breaking masses 
and usually exhibit certain hierarchy siructurebz3. In this work, in order to make our 
discussion as general as possible, we treat all ó;s as free parameter, and limit them by some 


physical observables. The constraints we consider include 


(I) The recently measured SM-like Higgs boson mass mp. In the MSSM, this mass is 
determined by the renormalized self-energies of the doublet CP-even Higgs fields, h, 
and hg, and the transition between them. Squarks contribute to these quantities 
through the g*qS and q*qSS interactions with S denoting either h, or dd . In 
the presence of the flavor mixings, both the interactions and the squark masses may 
be quite different from those in the case of 0;s = 0, and so is the SM-like Higgs boson 
mass. Among the flavor mixing parameters ó;, the Higgs mass is more sensitive to the 


chiral-flipping ones zpr and ógr. 


In this work, we get m, in the MSSM by the code FeynlliggsD3. In our scan over 
the parameter space of the low energy MSSM, we require the mass to be about 2GeV 
around its measured central value, i.e. 123GeV < mp € 127GeV. While for the MSSM 
in heavy SUSY case (see below), noting that the mass obtained by FeynHiggs suffers 
from potentially large theoretical uncertainties, we require a moderately wider range, 


ie. 121GeV < mp € 129GeV. 


(II) The LHC search for SUSY. By now both the ATLAS and CMS collaborations have 
paid great effects in searching for the signals of gluino, squarks as well as charginos 
and neutralinos, and based on certain assumptions, they exclude some SUSY particles 
up to about 1 TeV|28]. These obtained results, however, can not be applied directly 
to a general SUSY case, and in order to implement the LHC constraints, one has 
to perform detailed Monte Carlo simulation for each SUSY parameter point with the 
same strategies as those of the collaborations, then compare the simulated results with 
the LHC data|29]. In practice, such a process is rather time consuming, and can not 
be applied to an extensive scan over the SUSY parameter space, where a large number 


of samples are involved. 


In order to simplify our analysis, we note that by now gluino is preferred to be_at 
TeV scale without considering special cases such as compressed SUSY — 
while the second and third generation squark may still be as light as several hundred 
Gev ed lid. especially in the presence of the flavor mixing when the limitation on the 
squark spectrum can be further relaxed|23]. So we make following assumption in our 


discussion 
Mga Mu, 2 200GeV, mg >1TeV, mq, Mu, Mo, > 500GeV. (6) 
As will be shown below, our conclusions are not sensitive to such assumption. 


(IIT) The metastability of the vacuum state. This constraint reflects the fact that squarks 
as scalar fields contribute to SUSY potential, and consequently their soft breaking 
parameters should be limited by the stability (or more general metastability) of the 
vacuum state " 33. Assuming only zpr or dp, contributes to the potential, the 
metastability requireslai] 32 


Mê, + Mê, + M3 cos? B 


Mt 
fp) S 1.2 x ————————5———————, 
| HA V9 MP 
2 2 2 2 
EET m; Mg, + Mg, + M4 cos? B 
R| eo OTT 
V2 ME, 


|A| S 2.674/ M, + Mg, + M4 cos? B. (7) 


Note in previous study on top quark flavor changing neutral current process, this 


constraint was usually missed. 


(IV) EW precision observables My and sin? 6.7. In SUSY, the corrections to My and 
sin? 0,;; are dominated by squark loops, and their sizes reflect the mass disparity of 
left-handed SU(2) doublet squarks. To a good approximation, these two quantities 
are related with the óp parameter by 


Mw Gy . cys? 
Mw ~ — sp, sin? besp c -~“™ dp, 8 
wag ggi Osin? Bays = g bp (8) 


where 


j= a = (9) 


In this work, we repeat our previous calculation of 0M and 6 sin? bepp in Lo where 
three generation squarks are considered to implement the SU(2) relation between 


(M2)rr and (MẸ )LL; and require 
óMw <21MeV, — ósin?6,5; < 19.6 x 1077, (10) 


which are their allowed ranges at 2c level after considering experimental and theoret- 


ical uncertainties|33|. 


(V) Constraint from B — X,y. In the MSSM, the SUSY contributions to B + X,y come 
from four kinds of loops mediated by charged Higgs bosons, charginos, neutralinos 
and gluinos respectively. We calculate these contributions by the code FeynHiggs|27], 
and require 3.04 x 1074 € Br(B + X,y) € 4.02 x 1074 (corresponding its 2c allowed 


range by experiments|34]) in our parameter scan. 


Since the neutralino contribution is usually small, the expression of B — X,y in the 


NMSSM is roughly identical to that of the MSSM. 


About above constraints, it should be noted that constraints (I), (III) and (IV) do not 
diminish as SUSY scale becomes higher, that is, they are not decoupled in heavy SUSY 


limit; while the process B — X, does not possess such a property. 


FIG. 1: Feynman diagrams of the SUSY-QCD contribution to t > ch. If charm quark mass is 


neglected, the contribution from diagram (c) vanishes. 


III. SUSY PREDICTION ON THE RATE OF t ch 


In the MSSM, the dominant contribution to the process t — ch arises from the SUSY- 
QCD diagrams shown in Fig[l] The relevant Lagrangian is given by 


= vV2g,[9, ( -Uza Pr T Us, PR) dail i; jCj + g«(— Uza PL ae Usa Pr) doit izt] t h.c. 


6 
+ Y Copdah + Yitth, (11) 
a,B=1 


where g and ĝa denote gluino and squark in the mass eigenstate respectively, T is the 
Gell-Mann matrix with 7 and j representing color indices, U; is the 6 x 6 rotation matrix 


to diagonalize the mass matrix for up-type squarks, Cag parameterizes the coupling of the 


m cosa 
v sinf 


Higgs boson with squark mass eigenstates ĝa and qs, and Y, = with o being the 


rotation angle to diagonalize the CP-even Higgs mass matrix. The amplitude of t > ch can 


then be expressed by 
iM(t > ch) = ü(p;) [Fir + Fo) Pr, + (Fin + Fog) Pu]u(p;), (12) 


where, after neglecting charm quark mass, F;s are given by 


g? 


Fin — A [CagUs, Usem CoP? pi. Pt. T5. me, 


mi.) 
i og Uds UsgmCan (pz, Ph Dis M3, mi, m); 
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ur 
Fir = G3 Cas VeaU29mgCo(Ves Ph Pes m3; mz 49 à.) 


E TTT m mis ls 


dB 
igmg " 
Foy = tein YAUS, Uso Bo Pes m3, mi, ), 
d en ‘ 


In above expressions, p, p. and p, denote the momentums of top quark, charm quark, 
and Higgs boson respectively, mg and mg, represent the masses of gluino and squark re- 
spectively, and By , Co and Cj, are the standard two-point and three-point loop functions 
respectively|36]. In the heavy SUSY case discussed below, since the involved sparticle masses 
are much larger than mų, the contribution from Cj, can be safely ignored, and Bo, Co can 


be approximated by 


1 1 ô Oa 1 p? 
Cor i pt, mi, mi mày) = ar Ep np) - Fag meal + 005) 
g a a j 
à, 2 
Bo(pe,m3,mz,) = 1 1 n a + (55), (14) 
x: g 


where ĝa is defined as à, = m; /m;. 


A. t—> chin low energy SUSY 


As shown in bil. the SUSY-QCD contribution to t — ch in low energy MSSM has 


following features 


e In case that only one flavor mixing parameter ô; is non-zero, the rate of t — ch 
increases monotonously as the 6; becomes larger, while if several non-vanishing ó;s 


coexist, their effects may cancel each other out. 


e Since the effective hct vertex involves both chiral-flipping and flavor-changing, the 
chiral-conserving parameter ó7;/ógg must be accompanied with chiral-flipping ht7fg 
or hi5; interaction in contributing to the vertex, while the chiral-flipping parameter 
Órn/Ógr alone is able to lead into the hct vertex. As a result, the effective vertex is 
usually more sensitive to órg and dry, if we do not consider the constraints on the 


mixing parameters. 


e Unlike the other top quark FCNC processes, the Super-GIM mechanism does not apply 
to the decay t — ch. But since there exists a strong cancelation between diagram (a) 
and (b) in Fig[I] (see discussion below), the rate of t + ch depends on the soft mass 


parameters in a complex way. 


In the following, we do not intend to exhibit these features, but instead, noting that 


the constraints (I-III) introduced above were not considered before we try to figure out the 
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TABLE I: Benchmark points in low energy SUSY which correspond to very optimal cases in 
predicting the rate of t — ch. Points 1 and 2 are obtained from Scan-I, and Points 3 and 4 are 


from Scan-II. All these points satisfy the constraints listed in the text. 


EEUTRTETTPEFRETE 


order of Br(t — ch) that SUSY can predict after considering these constraints. For this 


purpose, we perform two independent scans over relevant SUSY parameters by imposing the 


constraints. Details of our scans are as follows: 


e Scan-I: We restrict our discussion in the MSSM, and calculate the Higgs mass with 


the code FeynHiggs. The parameter region we explore is given by 


500GeV < mg,,mqg,, mu,,mp,,mp, € 2TeV, 200GeV < my, < 2TeV, 
|A| € 6,/mg,;mu,, 1TeV € m; <2TeV, 1< tanp < 40, 
-lXÓór,ógg <1, —2.0 < dre, ÔrL € 20, —0.5 < d¢p, 5%, < 0.5, 


400GeV < my < 2TeV, —2TeV < u < 2TeV. (15) 
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In drawing up the strategy of this scan, we note that, although the down type squark 
parameters like 07, and Mp do not affect the rate of t — ch, they are needed in 
dp and B + X,y calculation. So to make our conclusions as general as possible we 
vary them in reasonable regions. We also note that since too many parameters are 
involved in the scan, the traditional random scan method is not efficient in searching 
for maximal value of Br(t — ch). So we adopt the Markov chain method in doing 
such a job. During the scan we adjust the optimal value by the results obtained from 


previous samplings until it reaches some stable values. 


e Scan-II: Same as Scan-I, but in order to relax the Higgs mass bound, we go beyond 
the MSSM by considering extra contribution to the mass. To be more specific, now we 
write the Higgs mass as m? = m? msgm + A? v? sin? 26 with A being a free parameter 
in the range from 0 to 0.7. This treatment of the Higgs mass is motivated by the 
singlet extensions of the MSSM such as the NMSSM, where the interaction between 
the singlet Higgs field and the doublet Higgs fields in the MSSM provides an additional 

Ri In order to get a large contribution from 


contribution to the mass at tree level 


this singlet extension, we set tan B = 1.5. 


In Table [I| we show two benchmark points for each scan. These points correspond to 
very optimal cases in predicting a large rate of t > ch. After analyzing our scan results, we 


have following observations 


e A small mg (around its experimental lower bound) seems to be favored to maximize 
the rate of t — ch after considering the constraints. Meanwhile, it is interesting to 
learn from Table [I] that, although we allow my, to be as low as 200GeV, the optimal 


points do not correspond to low Mys. 


e In low energy SUSY, the size of Br(t — ch) for the optimal points in Table [lis not 
sensitive to the parameters m, and u. For example, our results indicate that shifting 
ma from its value in Table[I] by 100 GeV only results in a change of Br(t > ch) by 
less than 196. While as will be shown below, these two parameters play an important 


role in determining the rate of the rare decay in heavy SUSY limit. 


e Among the considered constraints, the most stringent ones come from the Higgs mass 


and the LHC search for SUSY. As a result, the branching ratio of t — ch can only reach 
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1078 in low energy MSSM, which is about one order smaller than previous predictions 
obtained before the advent of the LEC. In scan-II, however, since the Higgs mass 


constraint on squark masses is comparatively relaxed, Br(t — ch) can reach 1075. 


About our results on t — ch, we remind that we do not include the SUSY-EW contri- 
bution to t — ch. The reason is the amplitude of the SUSY-EW contribution is roughly 


determined by omms« / Max(m;,,, m2.) from naive estimation B3, while the SUSY-QCD con- 


: 


tribution is determined by asmg/Max(mz,,, m; 


), where mg,, ma, and mz+ denote the mass 
scales for down-type squarks, up-type squarks and charginos respectively. Noting that m, 
is not much smaller than mg as suggested by the LHC search for the second and third gen- 
eration squarks bd bd, and also that the flavor mixings in the down-type squark sector are 
more tightly constrained by B-physics than those in up-type squark sector, we conclude that 
the SUSY-EW contribution should not be comparable with the SUSY-QCD contribution. 
So our estimates on the magnitude of t > ch will not change after including the SUSY-EW 
contribution. Another reason to neglect the EW contribution is, once considering it, too 
many parameters will be involved, but meanwhile this does not change our conclusion. 
Before we end this subsection, we have two comments. One is in extensions of the MSSM, 
the decay chain of a certain sparticle may be quite different from its MSSM prediction, and 
the analysis of the ATLAS and CMS collaborations in searching for SUSY may become 
irrelevant|38]. As a result, the constraint on the sparticle mass may be relaxed. This in 
return may push up the SUSY prediction on the rate of t — ch. Moreover, in extensions 
of the MSSM the couplings of the Higgs boson with quarks and squarks may be slightly 
changed. The influence of such changes on the rate of t — ch is usually not as signficant as 
the relax of the Higgs mass constraint. The other comment is that, given Br(t > ch) ~ 1075, 
it is difficult to detect such a top quark rare decay at the 14- TeV LHC with an integrated 
luminosity of 100 fo (Bo, but at future linear colliders such as TLEPug) detection of the 


decay is still possible. 


B. Remanent effect of t > ch in heavy SUSY limit 


As mentioned in Section I, because the hg*q¢ coupling strength is mainly determined by 
soft SUSY breaking parameters, the SUSY-QCD contribution to the effective Act interaction 


may exhibit remanent effect of SUSY in heavy SUSY limit. In order to investigate such an 
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FIG. 2: Dependence of Br(t — ch) on the common squark mass scale Msygsy. In getting this 
figure, we assume all soft masses equal to Msusy, ôLL = 0.7, tan 8 = 4, 10 and m4 = 400, 800GeV . 
We set A, = Msusy and u = 0 for the left panel, and A; = 0 and u = Msusy for the right panel. 


Note the lines corresponding to different choices of ma and tan f overlap in left panel. 


issue, we in the following assume a common SUSY mass scale mg, = ma, = mp, = Mu; 


mg = Msusy, and study the dependence of Br(t — ch) on Msusy for different choices of 
A, and u in Fig[2] In getting Fig[2] we consider the case that only 6,7 is non-vanishing and 
fix dz, = 0.7, tan B = 4,10 and m4 = 400, 800GeV. We set A; = Msusy and u = 0 for left 
panel, and A, = 0 and u = Msusy for right panel. These settings are only for exhibiting 
the decoupling behavior of Br(t — ch) and we do not consider the constraint from the 
LEP search for charginos, which requires u Z 103GeV. Fig[2| then indicates that SUSY has 
remanent effect on the rare decay rate only when p is at SUSY scale and meanwhile m4 is 
at weak scale, and the size of the effect depends strongly on m, and tan 6, e.g. small values 
of m4 and tan 8 tend to enhance the effect. We also investigate the case that only órg and 
dry are non-vanishing, and we do not find such remanent effect in heavy SUSY limit for any 
choices of u and mA. 

The behaviors shown in Fig[2] can be understood by the effective Lagrangian that de- 


scribes the Higgs and quark system. After including loop effects, one can write down the 
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Lagrangian as follows 


3 
L = Salmi; + ôm) Pray + ha (Y7 + SYL) Prd} + h«c. 


ij=1 


3 
3L Vim; + ôm) PL (Vra); + AVÀ): (Yy + 9Y5) PL(Vra);) + hc. 
ij=1 
: h : dm! 
= S{GmiPridi Y ,ümPay + 5 hq [Vg (0Y' — -p Velu Pog +h.c., (16) 


i=1 ij=l 
where q; and m;; = Yv are quark field and its mass matrix at tree level with 7, j denoting 
flavor indices, and 6m’ and Y” represent loop corrections to the mass matrix and the 
Yukawa coupling respectively. The second equation reflects the definition of quark mass 
eigenstate with Vz; and Vg denoting the rotation matrices for left-handed quarks and right- 
handed quarks respectively. After such a definition, the loop corrected mass matrix m’ + ôm’ 
is diagonal with its diagonal element m; representing physical quark mass determined by 
experiments. At this stage, the correction to the 4;q; interaction is given by Aq;|Vg(0Y' — 
óm'/v)Vr];;Prq; + h.c.. Obviously, if bY’ = óm'/v, new physics contribution to the Aq;q; 
interaction vanishes. In actual calculation, dm’ is obtained from q; — q; transition diagrams 
like diagram (b) of Fig] without the emission of the Higgs particle, and 6Y” comes from the 
vertex correction like diagram (a) of Fig[I] The effective Lagrangian then indicates that the 
two contributions should cancel out each other in contributing to the hqq, interaction. 

As far as the SUSY-QCD correction to the qq; interaction is concerned, its behavior 
in heavy SUSY limit can be analysised with the mass insertion approximation. In this 
method, squark masses are taken to be the diagonal elements of the squark mass matrix, 
and the non-diagonal elements are treated as interactions. In order to illustrate this method 
in explaining the remanent effect, we first consider the well studied SUSY-QCD correction 
to the Abb "T In this example, to get dm, one needs to insert the br, — bg transition 
by odd times into the sbottom propagator entered in bottom quark self-energy diagrams, 
and sum the corresponding contribution to infinite orders of the insertion. One can check 
that, with one more insertion, the corresponding contribution is suppressed by a factor 
mp(A, — utan 3)/M2ysy compared with that without the insertion, and only for the first 
order insertion, ôm, is not suppressed by Msysy, which means that the corresponding 
contribution is non-decoupled in heavy SUSY limit. In a similar way, one can check that 


even times of the insertion are needed to get the expression of óY, in calculating the Abb vertex 
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TABLE II: Same as Table[I] but for heavy SUSY case. Point 5 and Point 6 are taken from Scan-III 


and Scan-IV respectively. 


|. | Mga(GeV)| Mas(GeV) | Mya(GeV) | Mis (GeV) | A (GeV) "Em 
Point 5 4484 4039 6871 6839 6878  |-0.4885 1.202 -1.208 
Point 6 4158 4046 6886 6954 3944  —-0.3382| -0.1845 0.7258 


ÓRR mg(GeV) | ma(GeV) | u(GeV) Br(t > ch) 
Point 5| 0.7352 6.13 4095 800 -5943 EN 1.04 x 1078 
Point 6| 0.7734 4138 800 -18320 | 0.690 1.65 x 1078 


correction, and only for the zero-th insertion, the contribution is non-decoupled. Putting 


these two contributions together, one can learn that the non-decoupled terms proportional 
to A, is exactly canceled out, and the remaining contribution is proportional to the well 
known form pmg/Méysy (tan B + cota) ~ (—2mZums)/ (Màygy mA) tan B cos 28|[14]. This 
effect, albeit scaling as 1/m, does not diminish as jj ~ mg ~ Msusy approaches infinity 
for m4 at weak scale, and is therefore dubbed as the remanent effect of SUSY. 

Next we turn to analysis the SUSY-QCD contribution to the hct vertex. Since such 
a interaction involves both chiral-flipping and flavor-changing, appropriate insertions are 
needed to accomplish both the tasks. For the chiral-flipping mixings zr and Ógr, their 
role is quite similar to Aj in the SUSY-QCD correction to the hbb coupling, and their non- 
decoupling contribution is completely canceled out. While for the chiral-conserving mixings 
ÔLL and Opp, their contribution to the Abb vertex can be split into two parts with one part 
proportional to A; and the other part proportional to u. Fig[2] then reflects that the non- 
decoupling contribution of the former part is exactly canceled out, while that of the latter 
part is maintained if m4 is not at the same order as Msygy. This situation is actually similar 
to the SUSY-QCD correction to hbb vertex with the only difference coming from the fact 
that such remanent effect is not enhanced by tan 8. In fact, if both Msysy and mA approach 
infinity simultaneously, the genuine SUSY contribution to t — ch should vanish since now 
the Higgs sector of the MSSM is identical to that of the SM, while if only Msysy approaches 
infinity, the Higgs sector is described by a Two-Higgs-Doublet model, and SUSY may leave 
its imprint in Higgs sector]. Our results in FigD]actually reflect such a possibility. 

We finally discuss how large SUSY can predict the rate of t — ch in heavy SUSY case. For 
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this purpose, we require all squarks to be heavier than 3TeV and perform two independence 
scans over relevant SUSY parameter space by considering the constraints listed in Section 


II. These scans are 


e Scan-III: Similar to Scan-I except that we fix m4 = 800GeV and consider following 


parameter space 


ATeV < MQ, MQ, my,, mu,, MDa; mp, € TTeV, |A;| € 6,/mQ,mMu,, 
ATeV < m; € 10TeV, |u| € 20TeV, 1< tanp < 40, 


—l€ópnógg <1, —2.0 < ôrr, Ör, €20, —0.5 < 622,62; « 0.5. (17) 


In our scan, we do not consider very large squark soft breaking parameters because 
in such a case, the Higgs mass calculated by FeynHiggs suffers from large theoretical 


uncertainties. 
e Scan-IV: Similar to Scan-II except that the scan regions are now given by Eq.(17). 


For both scans, we find the rate of t + ch may reach 107? in optimal case, and a smaller 
value of m4 can lead to a larger branching ratio. In Table we provide two benchmark 
points for future study with Point 5 obtained from Scan-III, and Point 6 from Scan-IV. One 
can learn that, compared the prediction in heavy SUSY with that in low energy SUSY, 
although the optimal values of Br(t — ch) is suppressed by at least two orders, they are 
still 10° times larger than the corresponding SM prediction. 


IV. CONCLUSION 


In this work, we studied the top quark FCNC decay t — ch in the MSSM under the 
constraints from the Higgs mass measurement, the LHC searches for sparticles, the vacuum 
stability, the precision electro-weak observables and B — Xy. From a scan over the relevant 


parameter space, we found: 


e Due to the strong constraints from the measured Higgs mass and the results of SUSY 
searches at the LHC, the branching ratio of t > ch can only reach O(10 9) in the 
MSSM, which is about one order smaller than the old results. 
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e In the singlet extension of the MSSM, which can lift the Higgs mass at tree level, 


Br(t — ch) can reach O(10~°) in the allowed parameter space. 


e The chiral-conserving mixings ó;; and prr can induce SUSY remanent effect on the 
rate of t — ch. For heavy squarks and gluino above 3T'eV, Br(t — ch) can still reach 
1078. 
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